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ABSTRACT 

Hot gas ingestion problems for STOVL 
(Short Take-Off, Vertical Landing) aircraft are 
typically approached with empirical methods and 
experience. In this study, the hot gas environment 
around a STOVL aircraft was modeled as multiple 
jets in crossflow with inlet suction. The flow field 
was calculated with a Navier-Stokes, Reynolds- 
averaged, turbulent, 3D CFD code using a multi- 
grid technique. A simple model of a STOVL aircraft 
with four choked jets at 1000 K was studied at 
various heights, headwind speeds, and thrust splay 
angles in a modest parametric study. Scientific 
visualization of the computed flow field shows a pair 
of vortices in front of the inlet. This and other 
qualitative aspects of the flow field agree well with 
experimental data. 

NOMENCLATURE 

D, = characteristic length of jet nozzles, 0.0366 
m (1.44 in) 

H = distance from ground to bottom of the 
aircraft (aircraft altitude or height) 

U„ = headwind velocity 

V, = jet nozzle exit velocity, 633 m/s (2080 ft/sec) 

8 = thrust splay angle measured from the 

downward vertical toward symmetry plane 
k = turbulent kinetic energy 

£ = turbulent energy dissipation 

x = axial Cartesian coordinate, zero at upstream 

boundary 

x-y = vertical plane aligned in axial direction 
x-z = horizontal plane 

y = vertical Cartesian coordinate, zero at ground 
plane 

y-z = vertical plane aligned in spanwise direction 


z = spanwise Cartesian coordinate, zero at 

aircraft centerline plane 

INTRODUCTION 

Hot gas ingestion can cause significant 
problems for a STOVL (Short Take-Off, Vertical 
Landing) aircraft including reduced thrust and 
compressor stalls. These problems involve many 
hazards for the pilots including very hard landings. 
During the design of a STOVL aircraft, hot gas 
ingestion problems are typically approached with 
empirical methods and experience . Given the 
power of today's supercomputers and workstations, 
numerical methods employing efficient algorithms 
are becoming a viable engineering tool for analysis 
and design. In a previous endeavor, VanOverbeke 
& Holdeman 34 proved the feasibility of CFD analysis 
for hot gas ingestion. This study is a follow-on 
effort exploring the practicality of using an efficient 
numerical method for the problem of hot gas 
ingestion. 

FLOW FIELD DESCRIPTION 

Ingestion of hot gases generates problems 
in two ways: an average temperature rise results in 
a loss of engine thrust, and a temperature distortion 
may cause the engine to stall. Engine exhaust 
gases may be ingested by far-field and/or near-field 
mechanisms. A schematic of these mechanisms is 
shown in Fig. 1. 

The far-field mechanism results from the 
exhaust gases impinging on the ground and forming 
radial wall jets which flow forward, separate, and 
mix with the headwind. Near-field ingestion occurs 
with multiple jet configurations. Wall jets flowing out 
from the lift jets meet and create an upflow or 
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Hot gas ingestion mechanisms 

Figure 1 


fountain. This fountain flow can impinge on the 
aircraft’s underside, flow along the fuselage to the 
engine inlets, and be ingested. The gases ingested 
by this near-field mechanism tend to be hotter, 
giving greater temperature distortion than those 
ingested by the far-field mechanism. 

As stated earlier, the hot gas environment 
around a STOVL aircraft was modeled as multiple 
jets in crossflow with inlet suction. Mass sources 
represent the nozzle exits, and a mass sink at the 
end of the inlet provides the suction. The mass 
injected by the nozzles balances the mass removed 
by the inlet suction. This configuration derives from 
the previous study by VanOverbeke & Holdeman 34 . 
To meet the requirements of the CFD code, the 
aircraft model is placed in a confined flow, i.e., a 
'wind tunnel’. Also, the aircraft model has no angle 
of attack due to the use of a cartesian grid based 
flow solver. 

The STOVL aircraft model (see Fig. 2) is 
composed of rectangular solids for the fuselage and 
engine. For computational simplicity, the nose and 



the tail of the aircraft reach to infinity, and the 
model lacks wings. Baffles on the sides of the 
fuselage comprise the walls of the inlet. The 
nozzles are square in cross-section and are flush 
with the bottom of the aircraft. The square cross- 
section of the jets and the rectangular aircraft body 
result from the use of the cartesian grid. 

The four choked nozzles inject air at 1000 K 


(1340°F) straight down into the flowfield with a 
velocity of 633 m/s (2080 ft/sec). Each lift jet 
issues from the nozzle exit in a uniform flow. The 
headwind is also a uniform flow, but at a 
temperature of 300 K (81 °F). This approximates an 
aircraft landing with a forward speed or an aircraft 
facing into a wind which lacks a boundary layer. 
In the baseline case, the headwind (UJ flows at 3% 
of the jet velocity (V,), or about 19 m/s (37 kts), and 
the distance from the ground to the bottom of the 
aircraft model (H) is four times the characteristic 
length of the nozzles (4 Dj). The parametric studies 
include: various altitudes (H = 2 to 32 Dj) at a 
constant headwind (U„ = 0.03 V,); various headwind 
speeds (U„ = 0.01 to 0.09 V ( ) for a constant aircraft 
altitude (H = 4 Dj); and various thrust splay angles 
(8 = 0° to 45°) for a constant height (H = 4 D,) 
and headwind speed (U„ = 0.03 V,). 

The physical dimensions of the aircraft 
model are given in Table I. Note that the forward 
and aft nozzles have the same side-to-side 
separation, i.e., they are in-line, not offset. 


Table I 

STOVL aircraft model dimensions 


Fuselage: 

width (nose) 

s 

2.25 D. 

width (tail) 

= 

5.05 D, 

height 

=5 

2.5 D, 

length 

= 

oo 

Inlet: 

width 

S 

1.4 D, 

height 

= 

2.5 D, 

length 

a 

9.5 D. 

Jets separation: 
(center-to-center) 

side-to-side 

= 

3.25 D. 

fore & aft 

s 

6.0 D, 


NUMERICAL DESCRIPTION 
Calculation domain 

The grid geometry used for the baseline 
case is shown in Fig. 3. Exhibited are the 
centerline plane, the ground plane, and a vertical 
spanwise plane at the end of the domain as well as 
the aircraft model. The grid shows the high density 
of the calculation nodes in the region of the jets. 
For all calculations, symmetry assumptions allowed 
calculating only half of the physical domain. 

The other boundary conditions for the 
calculation domain include an inflow simulating the 
headwind for the domain face in front of the aircraft 
model and an outflow condition for the domain face 
behind the model. All flow properties are defined 
for the inflow condition. The outflow condition, in 
contrast, merely assumes the properties of the 
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Figure 3 



Figure 4 


axially nearest cells. The 
top, bottom, and remaining 
side of the domain are no- 
slip, stationary walls as are 
the aircraft surfaces. All the 
walls assume adiabatic 
conditions for the energy 
equation. The symmetry 
plane also has a symmetry 
condition for the energy 
equation. 

For the height 
variation, the grid contains 
211,200 cells arranged as 
follows: 100 cells in the x 
direction, 44 cells in the y 
direction, and 48 cells in the 
z direction. The physical 
dimensions of the baseline 
grid are about 135 D, long, 

29 D| high, and 40 Dj wide. 

The aircraft to ground 
distance was varied by 
elongating the cells 
underneath the aircraft. This 
facilitated the comparative 
analysis. 

The headwind speed 
variation used a slightly 
modified grid; 100 cells in x, 

44 cells in y, and 60 cells in 
z yielding 264,000 total cells. 

(See Fig. 4.) This grid has 
a greater length (177 D|) and 
greater width (59 D,) than the 
baseline grid. Also, the 
distance in front of the 
forward pair of jets is greater (152 Dj versus 76 D, 
for the baseline grid) to accommodate the long 
region of hot gas in front of the inlet in the U_ = 
0.01 V, case. No grid modification was needed to 
vary the headwind speed. 

Flow solver 

The flow field in this domain was calculated 
with a Navier-Stokes, Reynolds-averaged CFD code. 
References 5-7 describe this steady-state CFD code 
(and its techniques) for the three-dimensional 
analysis of turbulent elliptic flows in a Cartesian 
coordinate system. 

The CFD code (CART3D) solves the time- 
averaged Navier-Stokes or Reynolds equations. 
The k-e turbulence model provides closure. The 
governing equations include continuity, x-, y-, and z- 
momenta, energy, turbulent kinetic energy, and 
turbulent energy dissipation. These equations are 
solved using a block-implicit multigrid algorithm 
developed by Vanka. 

CART3D uses a hybrid differencing scheme 
on a staggered grid. This means that the code 


uses central differencing or upwind differencing 
depending upon the cell’s Reynolds number. Also, 
the scalar properties (density, pressure, etc.) are 
calculated at the cell volume centers while the 
velocities are solved at the centers of the cell faces. 

The multi-grid technique speeds 
convergence by solving the equations on sequential 
grids of different cell densities. The flow is 
initialized on the coarsest grid which gets refined by 
the multi-gridding. Dividing each cell on a grid into 
eight equal cells refines the grid for the next grid 
level. A V-cycle of sweeps on the various grid 
levels is performed until the solution converges on 
the finest grid. This technique speeds convergence 
by dampening out errors with the various levels of 
grid refinement. 

To determine convergence, the residuals are 
non-dimensionalized by an appropriate number, and 
then the maximum of all the residuals is compared 
to the tolerance criterion. The tolerance criterion 
used by this study is 1% for the finest grid. All test 
cases used the third grid level for the finest grid. 
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FLOW FIELD FEATURES 





x-z plane temperature contours 


H - 4 Dj, U, 


0.03 V, 


Figure 5 


Lift jets 


Lift jet particles 


temperature contours 


Select particle traces 
H = 4D r U„= 0.03 Vj 

Figure 6 


Fuselage 


l” <jt 


Forward vortex 


re contours on ground plane 


A short study of the 
features in the baseline case 
will help bring out the 
differences caused by varying 
the aircraft altitude. 

Fig. 5 displays the 
temperature contours in an x- 
z plane near the ground. 

These contours show the 
locations of the forward 
vortex pair and the two 
ground vortices generated by 
the interaction of the jets and 
the crossflow. The axis for 
the forward vortex pair is 
perpendicular to the plane 
shown in Fig. 5 while the 
axes for the ground vortices 
are parallel to the plane. 

The forward vortex pair is 
smeared by the steady-state 
calculations but still agrees 
well with the time-averaged 
experimental data 8 ' 10 . 

The particle traces in 
Fig. 6 reveal ingestion of 
exhaust gases. The particle 
traces from the jet region 
show the forward vortex and 
ingestion into the inlet. The 
particle traces starting at the 
inflow boundary show the 
headwind's deflection around 
the forward vortex. 

Temperature contours along 
the ground plane are also 
shown for clarity of position. 

Fig. 7 shows a three- 
dimensional temperature 
contour for the baseline 
case. This isotherm is for 
325 K (125°F), a reasonable 
upper limit on the 
temperature of the air 
reaching the engine. With 
the ambient flow at 300 K 
(81 °F), this represents a 
temperature rise of more 
than 25 K (45°F) for the fluid 
inside the isotherm. The 
inlet is almost completely 
obscured by the contour. 

Clearly, the engine is 
exposed to a considerable 
amount of hot gas from the 
engine exhaust. The bubble in front of the inlet 
reveals the location of the forward vortex. Note that 
the clipping of the isotherm in the right side of Fig. 


3D isotherm (T = 325 K) 

H = 4 Dj, U„ = 0.03 Vp T. = 300 K, T| = 1000 K 

Figure 7 

7 is due to wall effects which have no consequence 
on the hot gas ingestion. 
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AIRCRAFT ALTITUDE 
PARAMETER VARIATION 


The aircraft altitude varied from 2 Dj to 32 
Dj. The cases actually computed over this range (H 
= 2, 3, 4, 5, 6, 8, 12, 16, 24, and 32 Q) were 
carefully chosen to capture the changes in tne flow 
field features. Fig. 8 shows temperature contours 
for select aircraft altitudes (2 D,, 4 Q, 8 D r 12 Dj, 
and 32 D,) in two planes: the ground plane and a 
vertical plane passing through the jets near the jet 
centers. These cases show the changes in the flow 
field affecting hot gas ingestion over the range of 
variation. In each case, the fuselage is mostly 
hidden by the vertical plane of temperature 
contours. The major effects of aircraft altitude can 
be seen in this figure: the forward vortex changes 
in character, the amount of hot gas ingested is 
reduced, and the ground vortices decrease in size. 

In Fig. 9, the temperatures at the cells in 
front of the mass sink or 'engine face’ are plotted 
against the aircraft altitude. The temperatures 
shown are the minimum, the maximum, and a 
weighted average based on the cell volumes. The 
spread of the minimum temperature and the 
maximum temperature shows the temperature 
distortion at the engine face. This and the average 
temperature rise above ambient is plotted in Fig. 
10. At low altitudes, the distortion is obviously 



Temperature contours in select planes 
(U. = 0.03 V;) 

Figure 8 


U_ = 0.03 V ( 



Plot of engine face temperatures 
vs. aircraft altitude 

Figure 9 


extreme; and it quickly diminishes with increasing 
altitude. The average temperature shows a similar 
behavior. The non-monotonic behavior at H = 3 D, 
in all of these curves appears to be physical. Other 
calculations on different grids exhibit the same 
trends. 


U_ = 0.03 v ( 



Distortion and average rise 
vs. aircraft altitude 

Figure 10 


HEADWIND SPEED PARAMETER VARIATION 

The headwind speed was varied from 1 % to 
9% of the jet velocity (U_ = 0.01 to 0.09 VI in 1% 
increments for a constant altitude (H = 4 Dj). Fig. 
11 displays temperature contours in the same two 
planes as in Fig. 8: the ground plane and a vertical 
plane passing through the jets near the jet centers. 
Selected headwind speeds (1%, 3%, 5%, 7%, and 
9%) show the changes in the flow field for the 
varying headwind speed. Again, the fuselage is 
mostly hidden by the vertical plane of temperature 
contours. In Fig. 11 the calculated domain includes 
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Temperature contours in select planes 


(H = 4 D,) 

Figure 11 


unrealistically representing a windy vertical landing, 
but they might be relevant for a low speed landing. 
For choked jets, the 9% headwind represents about 
a 110 knot headwind which would either be a 
hurricane or a slower than normal landing speed for 
a conventional fighter aircraft. 



Distortion and average rise 
vs. headwind speed 

Figure 13 


a much larger region in front of the inlet in 
comparison with Fig. 8. This is due to a very long 
region of exhaust gas which extends in front of the 
inlet along the ground in the U„ = 0.01 V. case. 

The effect of headwind speed on tne engine 
face temperatures can be seen in Fig. 12. Note 
that the minimum temperature declines rather 
steadily with increasing headwind speed. The 
temperature distortion at the engine face varies 
weakly with headwind speed at the low speeds and 
is greatest for the 7% case as shown in Fig. 13. 
One should note that these high velocity ratios are 


H = 4 D| 



Headwind Speed - % of Vj 


Plot of engine face temperatures 
vs. headwind speed 

Figure 12 


THRUST SPLAY ANGLE PARAMETER VARIATION 

A technique used to help control hot gas 
ingestion is to splay the jets. By angling the lift 
jets, the relative strengths of the fountain, upwash, 
and vortices are changed thus changing the flow 
structures affecting ingestion. In this study, the 
splay angle (5) of the thrust is measured from the 
downward vertical inward to the centerline plane of 
the aircraft model. To vary the thrust splay angle, 
the component velocities on the jets changed to 
provide the required angle while keeping the speed 
of the jet constant. Thus the direction of the lift jets 
changed while the geometry of the aircraft model 
did not change. 

Splaying all jets 

For the first variation of the thrust splay 
angle, all four jets were splayed the same amount. 
The splay angle varied from 0° to 45° (8 = 0° to 
45°) in 5° increments for a constant height (H = 4 
Dj) and constant headwind speed (U„ = 0.03 Vj). 
Fig. 14 shows temperature contours in two planes 
similar to those displayed in Fig. 8 & Fig. 11: the 
ground plane and a vertical plane that is now on 
the inner side of the lift jets (instead of near the 
center). Selected splay angles (0°, 5°, 10°, 20°, 
25°, and 30°) show the changes in the flow 
structures due to varying the angle of all the jets. 
The most noticeable change occurs in the length of 
the hot gas region in front of the inlet. A less 
noticeable change is the increase of hot gas directly 
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Temp, contours - splay all jets 
(H = 4 D jf U„ = 0.03 Vj) 

Figure 14 


in front of the inlet for the 5° case. This actually 
causes an increase in hot gas ingestion over the 0° 
case, instead of reducing it. 

This and other effects of thrust splay angle 
for all the jets can be seen in Fig. 15 and Fig. 16. 
In Fig. 15, the average temperature first rises and 
then drops until a splay angle of about 20° where 
it starts rising again. A local maximum exists at the 
25° point before the average temperature flattens 


H = 4 D,, U. = 0.03 


min — f— max average 



K 250 


-5 0 5 10 15 20 25 30 35 40 45 50 

Splay Angle on All Jets - deg. 


Plot of engine face temperatures 
vs. splay angle on all jets 

Figure 15 


out at its lowest value Both the minimum and 
maximum temperatures follow the same behavior. 
The distortion and average temperature rise 
displayed in Fig. 16 show the same patterns. Note 
that for this configuration of model and altitude, the 
jets will converge at the ground for a thrust splay 
angle of 22°. 



Distortion and average rise 
vs. splay angle on all jets 

Figure 16 


Splaying forward jets only 

In the second variation, only the forward jets 
varied in thrust splay angle. The rear lift jets 

maintained a 0° splay angle. Again, the thrust 

splay angle varied from 0° to 45° (8 = 0° to 45°) in 
5° increments at the same altitude (H = 4 Dj) and 
headwind speed (U_ = 0.03 V) just as in the all jets 
splayed variation. Fig. 17 snows the temperature 
contours in the same two planes as in Fig. 14: the 
ground plane and a vertical plane on the inner side 
of the lift jets. Again, the selected splay angles (0°, 
5°, 10°, 20°, 25°, and 30°) show the changes in the 
flow structures due to varying the angle of the 
forward jets alone. The flow field changes are 
basically the same as in the cases with all the jets 
splayed. 

The ingestion effects of thrust splay angle 
for the forward jets can be seen in Fig. 18. The 
average temperature first rises and then drops until 
a splay angle of about 20° where it rises very 
slightly. A local maximum exists at 25° splay angle 
before the average temperature drops and flattens 
out at its lowest value. Both the minimum and 
maximum temperatures follow a similar behavior. 
The distortion and average temperature rise 
displayed in Fig. 19 show the same patterns, just 
as when all the jets are splayed. Overall, splaying 
the forward jets alone gives the same effects on hot 
gas ingestion as splaying all the jets. 
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Forward jets splayed 
(H = 4 D jt IL = 0.03 V,) 

Figure 17 


CONCLUSIONS/REMARKS 

In summary, the average of the engine face 
temperatures decreases with increasing height but 
is relatively unaffected by headwind speed. Engine 
face temperature distortion also decreases with 
increasing height, but unexpectedly increases with 
headwind speed until the forward vortex is behind 
the inlet face. The headwind speed variation 
reveals that (for a constant height) the engine face 
temperature is dominated by near field ingestion 
effects. 

As for the thrust splay angle variation, 
splaying the jets inward (for a constant height and 
headwind speed) first causes a rise and then a 
rapid decrease in hot gas ingestion, although a 
local maximum exists when the jets converge at the 
ground plane. Also, splaying the forward jets alone 
instead of all the jets gives almost the same 
reduction in hot gas ingestion without as large a 
thrust penalty. 

A comparison of the H = 4 Dj, U = 0.03 V, 
cases in the height parameter variation and the 
headwind speed parameter variation show the 
effects of the calculation domain on the flow field. 
The vertical walls of the narrower calculation 
domain definitely affect the ground vortices, but no 
differences exist in the temperatures reaching the 


H = 4 Dj, U_ - 0.03 ^ 



Splay Angle on Forward Jets - deg. 


Plot of engine face temperatures 
vs. splay angle on forward jets 

Figure 18 


Distortion and average rise 
vs. splay angle on forward jets 

Figure 19 


engine. Essentially, the tunnel walls in the 
calculation domain used for the height parameter 
variation are sufficiently far from the aircraft model 
so as to not affect the desired quantities (engine 
face temperatures). If the overall flow field is of 
primary interest, then the tunnel walls would have 
to be farther from the aircraft model. 

This study did not address the importance 
of the aircraft geometry (fuselage, wings, tails, etc.) 
in relation to the flow field. Only one aircraft model 
was used, and it was quite simplistic. 

The last conclusion from this study concerns 
the practicality of using an efficient CFD code for 
parameter variation studies. The turn-around time 
on a Cray-2 supercomputer and state-of-the-art 
workstations allows quick parameter changes. 
Typically, a Cray-2 supercomputer solved the flow 
field in about an hour with a turn-around time of a 
day. A dedicated IBM RS-6000 workstation can 
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solve the flow field in about 6 hours and can 
actually give shorter turn-around than the shared 
supercomputer. 
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